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ABSTRACT 



In this paper, we continue to build support for the proposal to use gamma-ray 
bursts (GRBs) as standard candles in constructing the Hubble Diagram at red- 
shifts beyond the current reach of Type la supernova observations. We confirm 
that correlations among certain spectral and lightcurve features can indeed be 
used as luminosity indicators, and demonstrate from the most up-to-date GRB 
sample appropriate for this work that the ACDM model optimized with these 
data is characterized by parameter values consistent with those in the concor- 
dance model. Specifically, we find that (fl m , ^a) ~ (0.30, 0.70), versus (0.27, 0.73) 
obtained from the 5-yr WMAP data. We also carry out a comparative analysis 
between ACDM and the i? h = ct Universe and show that the latter is a better fit 
to the GRB data. We find that the optimal ACDM model fits the GRB Hubble 
Diagram with a reduced xiof ~ 1-79, whereas the fit using = ct results in a 
Xdof ~ 1-66. In both cases, about 20% of the events lie at least 2a away from the 
best-fit curves, suggesting that either some contamination by non-standard GRB 
luminosities is unavoidable, or that the errors and intrinsic scatter associated 
with the data are being underestimated. 

Subject headings: cosmology: observations, redshift, theory; early universe; gamma-ray 
bursts: general 
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1. Introduction 

For a given class of sources whose luminosity is accurately known, one may construct 
a Hubble Diagram (HD) from the measurement of their distance versus redshift. Such a 
relationship can be a powerful tool for probing the cosmological expansion of the Universe, 
but only if these sources truly function as standard candles. The cosmic evolution depends 
critically on its constituents, so measuring distances over a broad range of redshifts can 
in principle place meaningful constraints on the assumed cosmology. The discovery of 
dark energy was made using this method, in which the sources — Type la supernovae — are 
transient, though with a well-defined luminosity versus color and light-curve shape 
relationships (Riess et al. 1998; Perlmutter et al. 1998, 1999; Garnavich et al. 1998; 
Schmidt et al. 1998). Of course, one must also assume that the power of distant explosions 
can be standardized against those seen at much lower redshifts. 

The use of Type la SNe has been quite impressive, so one may wonder why there 
would be a need to seek other kinds of standard candle. But the reality is that several 
important limitations mitigate the overall impact of supernova studies. For example, even 
excellent space-based platforms, such as SNAP (Scholl et al. 2004), cannot observe these 
events at redshifts > 1.8. And this is quite limiting because much of the most interesting 
evolution of the Universe occurred well before this epoch. In addition, the determination of 
the supernova luminosity cannot be carried out independently of the assumed cosmology, 
so the Type la SN data tend to be compliant to the adopted expansion scenario (Melia 
2012a). The fact that so-called "nuisance" parameters associated with the data need to be 
optimized along with the variables in the model itself weakens any comparative analysis 
between competing cosmologies. There is therefore much more to learn about the Universe's 
history than one can infer from Type la SNe alone. 

In recent years, several other classes of source have been proposed as possible standard 
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candles in their own right. Most recently, the discovery that high-z quasars appear to be 
accreting at close to their Eddington limit (see, e.g., Willott et al. 2010), has made it 
possible to begin using them to construct an HD at redshifts beyond ~ 6 (Melia 2012b). It 
has also been suggested that Gamma-ray Bursts (GRBs) may be suitable for constructing 
an HD at intermediate redshifts, 1 < z < 6, between the Type la SN and high-z quasar 
regions. 

The possible use of GRBs as standard candles started to become reality after Ghirlanda 
et al. (2004) found a rather tight correlation between the peak of the gamma-ray spectrum 
E p (in a power vF v versus frequency plot) and the collimation-corrected energy emitted 
in gamma-rays. Since E 1 = _E 7iiso (l — cos#), where 9 is the jet half-opening angle, one may 
reliably estimate the isotropically equivalent energy -E 7i i SO and use this to infer a distance. 

Earlier, Schaefer (2003) had constructed the first GRB HD based on nine events using 
two luminosity indicators, and this was followed by Bloom et al. (2003a), who published 
a GRB HD with 16 bursts, assuming that the burst energy is a constant after correcting 
for the beam angle. These and other attempts at constructing a GRB HD were made with 
only a small fraction of the available data, using only one or two luminosity indicators. 
Unfortunately, all of them had error bars that were too large to provide useful constraints 
on cosmology. 

The feasibility of using this method became better grounded following the Ghirlanda 
et al. (2004a) discovery, when Dai et al. (2004) found tight constraints on the cosmological 
parameters using their correlation. Unfortunately, because of the current poor information 
on low-z GRBs, the Ghirlanda relation necessarily depends on the assumed cosmology. 
However, Dai et al. (2004) assumed a universal correlation independent of the expansion 
scenario, so their results could not be used for comparative studies. Following this early 
work, other authors attempted to circumvent the circularity problem by using a correct 
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approach (Ghirlanda et al. 2004b; Firmani et al. 2005; Xu et al. 2005; Ghirlanda et 
al. 2006; Wang et al. 2011). In the end, however, as is also true for Type la SNe, the 
correlation ought to be recalibrated for each different model because the best-fit correlation 
depends on the cosmology adopted to derive the burst luminosities. 

Of course, even with the emergence of more precise luminosity indicators, one must 
still deal with several significant challenges when trying to use GRBs to construct an HD. 
The luminosity of these bursts, calculated assuming isotropy, spans about 4 orders of 
magnitude (Frail et al. 2001). However, there is strong observational evidence (e.g., the 
achromatic break in the afterglow lightcurve) that the burst emission is collimated into a 
jet with the aforementioned aperture angle 9 (Levinson & Eichler 1993; Rhoads 1997; Sari 
et al. 1999; Fruchter et al. 1999). When one corrects for the collimation factor (1 — cos#), 
the gamma-ray energy tends to cluster around E 1 ~ 10 51 ergs, but the dispersion (~ 0.5 
dex) is still too large for these measurements to be used for cosmological purposes. This is 
why much effort has been expended since 2004 in finding other indicators from the GRB 
spectrum that provide more precise constraints on the luminosity. 

One useful application of these ideas involves the use of E p and t h (the so-called jet 
break time) as a measure of to determine -E 7i i S o (Liang & Zhang 2005). In this paper, we 
will follow this approach with two distinct goals in mind. First, several new GRB events 
have been detected in recent years that have spectral and lightcurve features (such as t^) 
with sufficient quality to help improve the previously assembled correlations. Second, and 
foremost, we wish to use this relatively new probe of the Universe's expansion to directly 
test the R h = ct Universe (Melia 2007; Melia & Shevchuk 2012) against the data and to see 
how its predictions compare with those of the ACDM cosmology. 

In the next section, we will describe the data we will use and our method of analysis. 
We will then first assemble the GRB HD in the context of the standard model, ACDM, 
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and demonstrate that the best-fit parameters obtained by fitting its predicted luminosity 
distance to the GRB observations very closely mirror those obtained through the analysis 
of Type la SNe. In § 4, we will introduce the = ct Universe and provide a brief overview 
of its current status as a viable cosmology. We will then construct the GRB HD for this 
expansion scenario, which is much easier to do than for ACDM because the former has only 
one free parameter — the Hubble constant H . Finally, we will directly compare the results 
of our fits to the data with both ACDM and i?h = ct. 



2. Observational data and Methodology 

Our GRB sample includes 28 bursts with a measurement of the redshift z, the spectral 
peak energy E p , and the jet break time % seen in the optical afterglow. All of these data 
have been obtained from previously published studies. Our complete sample is shown in 
Table 1, which includes the following information for each GRB: (1) its name; (2) the 
redshift; and various spectral fitting parameters, including (3) the spectral peak energy E p 
(with corresponding error &e p ), (4) the low-energy photon index a, (5) the high-energy 
photon index j3; (6) the 7-ray fluence S 1 (with error o"s 7 ); (7) the observed energy band; 
and (8) the jet break time % (with error a tb ). 

With the data listed in Table 1, we calculate the isotropic equivalent gamma-ray energy 
(-E 7l iso) using 

AnDl(z)S 

£ 7 ,iso - {1 + z) K, (1) 

where is the measured gamma-ray fluence, D L (z) is the luminosity distance at redshift 
z, and K is the i^-correction factor used to correct the gamma-ray fluence measured 
within the observed bandpass (taken to be 1 — 10 4 keV in this paper) and shift it into the 
corresponding bandpass seen in the cosmological rest frame. 



Table 1. GRB Prompt Emission Parameters and the Jet Break Time 



GRB 


z 


E P (*E P ) 


a 






Band 


*b(ot b ) 


References 






(keV) 






(10~ 6 ergs cm~ 2 ) 


(keV) 


(days) 




970828 


0.9578 


297.7 ± 59.5 


-0.70 


-2.07 


96 ± 9.6 


20 


- 2000 


2.2 ± 0.4 


1,2, 2,2 


980703 


0.966 


254 ± 50.8 


-1.31 


-2.40 


22.6 ± 2.3 


20 


- 2000 


3.4 ± 0.5 


3, 4, 4, 5 


990123 


1.6 


780.8 ± 61.9 


-0.89 


-2.45 


300 ± 40 


40 


- 700 


2.04 ± 0.46 


6, 7, 7, 6 


990510 


1.62 


161.5 ± 16.1 


-1.23 


-2.70 


19 ± 2 


40 


- 700 


1.6 ± 0.2 


8, 7, 7, 9 


990705 


0.8424 


188.8 ± 15.2 


-1.05 


-2.20 


75 ± 8 


40 


- 700 


1 ± 0.2 


» » . 2 


990712 


0.43 


65 ± 11 


-1.88 


-2.48 


6.5 ± 0.3 


40 


- 700 


1.6 ± 0.2 


8, 7, 7, 10 


991216 


1.02 


317.3 ± 63.4 


-1.23 


-2.18 


194 ± 19 


20 


- 2000 


1.2 ± 0.4 


11,4, 4, 12 


000926 


2.07 


100 ± 7 


-1.10 


-2.43 


26 ± 4 


20 


- 2000 


1.74 ± 0.11 


13, 14, 13, 15 


010222 


1.48 


291 ± 43 


-1.05 


-2.14 


88.6 ± 1.3 


40 


- 700 


0.93 ± 0.15 


16, 17, 17, 16 


011211 


2.14 


59.2 ± 7.6 


-0.84 


-2.30 


5 ± 0.5 


40 


- 700 


1.56 ± 0.02 


18, 19, 18, 20 


020124 


3.2 


86.9 ± 15 


-0.79 


-2.30 


8.1 ± 0.8 


2 - 


400 


3 ± 0.4 


21, 22, 22, 23 


020405 


0.69 


192.5 ± 53.8 


0.00 


-1.87 


74 ± 0.7 


15 


- 2000 


1.67 ± 0.52 


24, 24, 24, 24 


020813 


1.25 


142 ± 13 


-0.94 


-1.57 


97.9 ± 10 


2 - 


400 


0.43 ± 0.06 


25, 22, 22, 25 


021004 


2.332 


79.8 ± 30 


-1.01 


-2.30 


2.6 ± 0.6 


2 - 


400 


4.74 ± 0.14 


26, 22, 22, 27 


021211 


1.006 


46.8 ± 5.5 


-0.86 


-2.18 


3.5 ± 0.1 


2 - 


400 


1.4 ± 0.5 


28, 22, 22, 29 


030226 


1.986 


97 ± 20 


-0.89 


-2.30 


5.61 ± 0.65 


2 - 


400 


1.04 ± 0.12 


30, 22, 22, 31 


030328 


1.52 


126.3 ± 13.5 


-1.14 


-2.09 


37 ± 1.4 


2 - 


400 


0.8 ± 0.1 


32, 22, 22, 33 


030329 


0.1685 


67.9 ± 2.2 


-1.26 


-2.28 


163 ± 10 


2 - 


400 


0.5 ± 0.1 


34, 22, 22, 35 


030429 


2.6564 


35 ± 9 


-1.12 


-2.30 


0.85 ± 0.14 


2 - 


400 


1.77 ± 1 


36, 22, 22, 37 


041006 


0.716 


63.4 ± 12.7 


-1.37 


-2.30 


19.9 ± 1.99 


25 


- 100 


0.16 ± 0.04 


2, 2, 2, 38 


050401 


2.9 


128 ± 30 


-1.00 


-2.45 


19.3 ± 0.4 


20 


- 2000 


1.5 ± 0.5 


39, 39, 39, 40 


050408 


1.2357 


19.93 ± 4 


-1.98 


-2.30 


1.9 ± 0.19 


30 


- 400 


0.28 ± 0.17 


2, 2, 2, 41 


050416A 


0.653 


17 ± 5 


-1.01 


-3.40 


0.35 ± 0.03 


15 


- 150 


1 ± 0.7 


39, 39, 39, 40 


050525A 


0.606 


79 ± 3.3 


-0.99 


-8.84 


20.1 ± 0.5 


15 


- 350 


0.28 ± 0.12 


42, 42, 42, 42 
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Table 1 — Continued 



GRB 


z 


E P (a Ep ) 


a 






Band 


*b(ot b ) 


References 






(keV) 






(10~ 6 ergs cm -2 ) 


(keV) 


(days) 




060206 


4.048 


75.5 ± 19.4 


-1.06 




0.84 ± 0.04 


15 - 150 


2.3 ± 0.11 


39, 39, 39, 40 


060526 


3.21 


25 ± 5 


-1.10 


-2.20 


0.49 ± 0.06 


15 - 150 


2.77 ± 0.3 


39, 39, 39, 43 


060614 


0.125 


49 ± 40 


-1.00 




22 ± 2.2 


15 - 150 


1.38 ± 0.04 


39, 39, 39, 44 


070125 


1.547 


367 ± 58 


-1.10 


-2.08 


174 ± 17 


20 - 10000 


3.8 ± 0.4 


39, 39, 39, 45 



References-. rp^ c re f erences a pp ear j n the following order: redshift, spectral parameters, fluence, and break time: 
(1) Djorgovski et al. (2001); (2) Wang & Dai et al. (2006); (3) Djorgovski et al. (1998); (4) Jimenez et al. (2001); 
(5) Frail et al. (2003); (6) Kulkarni et al. (1999); (7) Amati et al. (2002); (8) Vreeswijk et al. (2001); (9) Stanek 
et al. (1999); (10) Bjornsson et al. (2001); (11) Djorgovski ct al. (1999); (12) Halpcrn et al. (2000); (13) Amati 
et al. (2006); (14) Xiao et al. (2009); (15) Sagar et al. (2001); (16) Galama et al. (2003); (17) Guidorzi et al. 
(2011); (18) Holland et al. (2002); (19) Amati (2003); (20) Jakobsson et al. (2003); (21) Hjorth ct al. (2003); (22) 
Sakamoto et al. (2005); (23) Berger ct al. (2002); (24) Price et al. (2003); (25) Barth ct al. (2003); (26) Moller et 
al. (2002); (27) Holland et al. (2003); (28) Vreeswijk et al. (2003); (29) Holland et al. (2004); (30) Grciner et al. 
(2003); (31) Klose et al. (2004); (32) Martini et al. (2003); (33) Andersen et al. (2003); (34) Bloom ct al. (2003b); 
(35) Berger et al. (2003); (36) Weidinger et al. (2003); (37) Jakobsson et al. (2004); (38) Stanck ct al. (2005); (39) 
Ghirlanda et al. (2008); (40) Ghirlanda ct al. (2007); (41) Godct ct al. (2005); (42) Blustin ct al. (2006); (43) Dai 
et al. (2007); (44) Delia Valle et al. (2006); (45) Chandra ct al. (2008). 
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Both ACDM and R h = ct are Friedmann-Robertson- Walker (FRW) cosmologies, but 
the former assumes specific constitutents in the density, written as p = p r + p m + p A , where 
p r , p m and pa are, respectively, the energy densities for radiation, matter (both luminous 
and dark) and the cosmological constant. These densities are often written in terms of 
today's critical density, p c = 3c 2 Hq/8tiG, represented as Q m = p m j p c , Q r = p r /p c , and 
Q A = pa/p c . In a flat universe with zero spatial curvature, the total scaled energy density 
is Q = Q m + Q r + = 1- In Rh — ct, on the other hand, the only constraint is the total 
equation of state p = wp, where w = —1/3. Later in this paper, we will discuss how these 
two formulations are related to each other, particularly the fact that w = —1/3 uniquely 
forces Vl m = 0.27 in ACDM when p A = -p A (Melia 2012c). 

In ACDM, the luminosity distance is given as 



Dt CDM (z) = f ^= sinn 

H VI | 



L a/2 x r dz \ , 2) 

\ Jo ^{i + z) 2 (i + n m z)-z(2 + z)n A / ' 

where c is the speed of light, and H is the Hubble constant at the present time. In 
this equation, Q k is defined similarly to Vt m and represents the spatial curvature of the 
Universe — appearing as a term proportional to the spatial curvature constant k in the 
Friedmann equation. Also, sinn is sinh when Qk > and sin when < 0. For a flat 
Universe with £l k = 0, Equation (2) simplifies to the form (1 + z)c/H times the integral. 
For the i?h = ct Universe, the luminosity distance is given by the much simpler expression 

D^= ct = ^-(l + z)Ml + z). (3) 

The factor c/Hq is in fact the gravitational horizon Rh(to) at the present time, so we may 
also write the luminosity distance as 

D^= ct = R h (t )(l + z)\n(l + z). (4) 

For each of the GRB sources listed in Table 1, we have derived the equivalent isotropic 
energy according to Equation (1) and listed it in Table 2. Here, E^^ M and E^\~ ct are the 
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Table 2. Derived Rest-frame Burst Properties in ACDM and R h = ct 



GRB \ogE> p (a E ,) log i^) log E$$ M (* E J log E*£*{v E J 
(keV) (day) (ergs) (ergs) 



970828 


2 


.77 


± 





.09 


0.051 ± 0.079 


53 


,48 


± 





,04 


53 


.38 


± 





.04 


980703 


2 


,70 


± 





.09 


0.238 ± 0.064 


52. 


,86 


± 





.04 


52 


.75 


± 





.04 


990123 


3 


.31 


± 





.03 


-0.105 ± 0.098 


54. 


,61 


± 





.06 


54 


.50 


± 





.06 


990510 


2 


.63 


± 





.04 


-0.214 ± 0.054 


53. 


.31 


± 





.05 


53 


.20 


± 





.05 


990705 


2 


.54 


± 





.03 


-0.265 ± 0.087 


53. 


.41 


± 





.05 


53 


.31 


± 





.05 


990712 


1 


.97 


± 





.07 


0.049 ± 0.054 


51 


.92 


± 





.02 


51 


.84 


± 





,02 


991216 


2 


.81 


± 





.09 


-0.226 ± 0.145 


53. 


,84 


± 





.04 


53 


.74 


± 





,04 


000926 


2 


.49 


± 





.03 


-0.247 ± 0.027 


53. 


.54 


± 





.07 


53 


.43 


± 


0. 


.07 


010222 


2 


.86 


± 





.06 


-0.426 ± 0.070 


53. 


.97 


± 





.01 


53 


.85 


± 





.01 


011211 


2 


.27 


± 





.06 


-0.304 ± 0.006 


53. 


.02 


± 





.04 


52 


.91 


± 


0. 


.04 


020124 


2 


.56 


± 





.07 


-0.146 ± 0.058 


53. 


.40 


± 





.04 


53 


.31 


± 


0. 


.04 


020405 


2 


.51 


± 





.12 


-0.005 ± 0.135 


53. 


.13 


± 





.00 


53 


.03 


± 


0. 


.00 


020813 


2 


.50 


± 





.04 


-0.719 ± 0.061 


54. 


.14 


± 





.04 


54 


.03 


± 





.04 


021004 


2 


.42 


± 





.16 


0.153 ± 0.013 


52. 


.67 


± 





.10 


52 


.57 


± 





.10 


021211 


1 


.97 


± 





.05 


-0.156 ± 0.155 


52. 


,14 


± 





.01 


52 


.03 


± 





.01 


030226 


2 


.46 


± 





.09 


-0.458 ± 0.050 


52. 


,90 


± 





.05 


52 


.79 


± 





.05 


030328 


2 


.50 


± 





.05 


-0.498 ± 0.054 


53. 


,58 


± 





.02 


53 


.46 


± 





.02 


030329 


1 


.90 


± 





.01 


-0.369 ± 0.087 


52. 


.20 


± 





.03 


52 


.16 


± 





.03 


030429 


2 


.11 


± 





.11 


-0.315 ± 0.245 


52. 


.26 


± 





.07 


52 


.16 


± 





.07 


041006 


2 


.04 


± 





.09 


-1.030 ± 0.109 


53. 


.01 


± 





.04 


52 


.91 


± 





.04 


050401 


2 


.70 


± 





.10 


-0.415 ± 0.145 


53. 


,64 


± 





.01 


53 


.55 


± 





.01 


050408 


1 


.65 


± 





.09 


-0.902 ± 0.264 


52. 


,41 


± 





.04 


52 


.29 


± 





.04 


050416A 


1 


.45 


± 





.13 


-0.218 ± 0.304 


50. 


,96 


± 





.04 


50 


.86 


± 


0. 


.04 


050525A 


2 


.10 


± 





.02 


-0.759 ± 0.186 


52 


,38 


± 





.01 


52 


.29 


± 


0. 


.01 
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Table 2 — Continued 



GRB 


log E' p (a E 


log t' h (a t , b ) 




i°g£ 7 fl L7V^) 




(keV) 


(day) 


(ergs) 


(ergs) 


060206 


2.58 ± 0.11 


-0.341 ± 0.021 


52.77 ± 0.02 


52.71 ± 0.02 


060526 


2.02 ± 0.09 


-0.182 ± 0.047 


52.43 ± 0.05 


52.35 ± 0.05 


060614 


1.74 ± 0.35 


0.089 ± 0.013 


51.25 ± 0.04 


51.22 ± 0.04 


070125 


2.97 ± 0.07 


0.174 ± 0.046 


53.99 ± 0.04 


53.87 ± 0.04 
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isotropic energies in ACDM and Rh = ct, respectively, remembering that all of the standard 
candle features must be re-calibrated for each assumed expansion scenario. Note, however, 
that E' p and t' h depend only on redshift z and are therefore independent of the assumed 
cosmology. One notices immediately that, although the numbers differ slightly, they are in 
fact remarkably similar, even though £)A CDM ari d D^ h ~ ct have quite different formulations. 
This is another consequence of the fact that ACDM mimics the = ct Universe quite 
closely, as we will discuss later in this paper (see also Melia 2012c). 

As we alluded to in the introduction, our approach is similar to that presented by 
Liang & Zhang (2005), in which we seek an empirical relationship between E 1;iso , E' , and 
t' h , known as the Liang-Zhang relation. Our form of the luminosity correlation is written as 
follows: 

log £ 7jiso = k + ki log E' p + k 2 log t' h , (5) 

where E' p = E p (l + z) in keV and t' h — t h /(l + z) in days. To find the best-fit coefficients k , 
Ki and k 2 , we follow the technique described in D'Agostini (2005). Let us first simplify the 
notation by writing x\ = \ogE' p , x 2 = logi b , and y = \ogE ltiso . Then, the joint likelihood 
function for the coefficients k , ki, k 2 and the intrinsic scatter <7 int , is 



L(k , «i, « 2 , O-int) OC Yl / 1 X 



cxp 



(6) 



2(^ + + 
where % is the corresponding serial number of each GRB in our sample. 

The best-fit luminosity correlation is shown in Figure 1, together with the data, for 
both ACDM (left panel) and R h = ct (right panel). For this exercise, we assumed a flat 
ACDM cosmology with Q m = 0.27 and Hq = 70 km s _1 Mpc -1 , obtained from the 5-yr 
WMAP data (Komatsu et al. 2009). Using the above optimization method, we find that in 
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ACDM the best-fit correlation between i? 7i i so , and E' and t' h , is 



log£ 7 , iso = (48.17 ±0.64) + (1.93 ±0.25) log E' p - (0.91 ± 0.34) log t' h 



(7) 



with an intrinsic scatter cr int = 0.27 ± 0.08. (In Table 2, this energy is labeled E^® M .) The 
best-fitting curve is plotted in Figure la. 
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Fig. 1. — The E 1 j ao versus E' p — t' h correlation in ACDM and = ct. The solid curves show 
the best fitting results from Equation (5). 

In the i?h = ct Universe, there is only one free parameter — the Hubble constant H . We 
note, however, that both the data and the theoretical curves depend on 1/H , since we do 
not know the absolute value of the GRB luminosity. As such, though formally Ho is a free 
parameter for both ACDM and the i?h = ct Universe, in reality the fits we discuss in this 
paper do not depend on its actual value. For the sake of consistency, we will adopt the 
standard H = 70 km s _1 Mpc -1 throughout our analysis and discussion. 



Using the above methodology, we find that the best-fitting correlation between E. 
and E' p and t' h , is now 



7, iso) 



log £ 7>iso = (48.13 ± 0.63) + (1.91 ± 0.24) log E' - (0.89 ± 0.32) \ogt' h 



(8) 
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with an intrinsic scatter a- m t = 0.26 ± 0.08. The best-fitting curve is plotted in Figure lb 
(and is labeled E^~ ct in Table 2). These coefficients are quiet similar to those obtained for 
ACDM. 

3. Optimization of the Model Parameters in ACDM 

The dispersion of the empirical relation for £7 7j is is so small that it has served well as 
a luminosity indicator for cosmology (Liang & Zhang 2005; Wang & Dai 2006). However, 
since this luminosity indicator is cosmology- dependent, we cannot use it to constrain the 
cosmological parameters directly. In order to avoid circularity issues, we use the following 
two methods to circumvent this problem: 




Fig. 2. — Plot of — log(likelihood) (left panel) and the intrinsic scatter a- m t (right panel), as 
functions of fl m , obtained by fitting the correlation with the joint likelihood method in a flat 
universe. 

Method I. We repeat the above analysis while varying the cosmological parameter Q m , 
though first under the assumption that the Universe is flat (Amati et al. 2008; Ghirlanda 
2009). The two panels in Figure 2 show that the values of — log(likelihood) and the intrinsic 
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scatter a- mt are indeed sensitive to Q m , showing a clear minimum around Vt m ~ 0.30. 
Moreover, the correlation slopes Ki and k 2 are also sensitive to the assumed cosmology, 
as shown by the two panels in Figure 3. Using the probability density function, the joint 
likelihood method allows us to constrain Q m to lie within the range 0.22 — 0.80 at the la 
confidence level. 




Fig. 3. — The slopes Ki (left panel) and k 2 (right panel) as functions of Q m , obtained by 
fitting the correlation with the joint likelihood method in a flat universe. 

If we release the flat universe constraint and allow Q m and Qa to vary independently 
(see Figure 4), the contours show that fl m and are poorly constrained; only an upper 
limit of ~ 0.9 can be set at lcr for both Q m and ^a- However, if we consider only a flat 
universe, the allowed region at the la level is restricted by the flat Universe (dashed) line 
and the la contour, for which 0.14 < Q m < 0.79 and 0.21 < Q\ < 0.86. The most probable 
values of Q m and I^a are (0.32,0.68). 

Method II. The distance modulus of a GRB is defined as 



^ = 51ogpx/10j?c) , (9) 
in terms of the luminosity distance D^. Using the Liang-Zhang relation, we can recast this 
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Fig. 4. — Contour confidence levels of Vt m and Qa, obtained by fitting the correlation with 
Method I. 

in the form 



However, since the luminosity correlation is cosmology-dependent, (x also depends on the 
adopted expansion scenario. We use the following approach to circumvent this difficulty 
(see also Liang & Zhang 2005). This procedure is based on the calculation of the probability 
function for a given set of cosmological parameters (denoted by Q, which includes both fl m 
and ^a) : 

Step 1. For a given cosmological model, we calibrate and weight the luminosity 
indicator corresponding to each choice of parameters Q. In each case, we calculate the 
correlation E 7 ^ so (fl; E' ,t' h ), and evaluate the probability [if(fj)] of this relation being the 
optimal cosmology-independent luminosity indicator via \ 2 statistics, i.e., 



p, = 2.5[k + «i log 4 + k 2 \ogt' h - \og(AnS y K) + log(l + z)\ - 97.45 . 



(10) 




N 



(11) 



The probability is then 




(12) 
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Step 2. We regard the correlation derived for each set of parameters as a cosmology- 
independent luminosity indicator without considering its systematic error, and calculate the 
distance modulus fiifl) and its error c^, given by 

1/2 



2.5 
In 10 



+ 



OS. 



7> l 



Ka 



+ 



1 + Zi 



(13) 



Since both (aKjKi) 2 and [c Zi /(l + Zi)] 2 are significantly smaller than the other terms in 
Equation (13), we ignore them in our calculations. 



a 




0.8 1.0 



Fig. 5. — Contour confidence levels of fl m and inferred from the current GRB sample, 
using method II. 



Step 3. We calculate the theoretical distance modulus /i(f2) for a set of cosmo logical 
parameters (denoted by ft), and then obtain x 2 from a comparison of with //(O), i.e., 



(14) 



Step 4- We then calculate the probability that the cosmological parameter set Q is the 
correct one according to the luminosity indicator derived from the cosmological parameter 
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set f2, i.e., we calculate 

p(ti | Q) oc e-* 2 ^' 2 . (15) 

Step 5. Finally, we integrate Q over the full cosmological parameter space to get the 
final normalized probability that the cosmological parameter set Q is the correct one, i.e., 

, x Lw(n) P (n I n) dd 

Jn w ( Q ) dQ 

Figure 5 shows the la to 3a contours of the probability in the (fi m , fi A ) plane. The 
contours show that at the la level, 0.13 < Q m < 0.51, but f2 A is poorly constrained; only an 
upper limit of ~ 0.82 can be set at this confidence level. However, if we consider only a flat 
Universe, then the allowed range of parameter space is limited by the flat Universe dashed 
line and the la contour, for which 0.26 < Q m < 0.49 and 0.51 < Q A < 0.73. The best fit 
values are (ft,fi A ) = (0.30,0.70). 



4. The R h = ct Universe 

In the previous section, we considered how the currently available sample of GRB 
events with spectral and lightcurve characteristics appropriate for cosmological work may 
be used to constrain the principal parameters of the standard model. The = ct Universe, 
on the other hand, has only one free parameter — the Hubble constant Hq. However, as we 
have already noted, none of the results presented in this paper depend on this constant, 
since 1/H enters into the determination of both the data and the theoretical curves. There 
is therefore no need to reproduce the kind of parameter optimization for = ct that was 
carried out for ACDM in § 3. But before we proceed to compare the Hubble diagrams for 
the i?h = ct Universe and the optimized ACDM model, we will first briefly summarize the 
Rh = ct cosmology, which is not yet as well known as ACDM. 
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One may look at the expansion of the Universe in several ways. From the perspective 
of the standard model, one guesses the constituents and their equation of state and then 
solves the dynamical equations to determine the expansion rate as a function of time. The 
second is to use symmetry arguments and our knowledge of the properties of a gravitational 
horizon in general relativity (GR) to determine the spacetime curvature, and thereby the 
expansion rate, strictly from just the value of the total energy density p and the implied 
geometry, without necessarily having to worry about the specifics of the constituents that 
make up the density itself. This is the approach adopted by i?h = ct. In other words, 
what matters is p and the overall equation of state p = wp, in terms of the total pressure 
p and total energy density p. In ACDM, one assumes p = p m + p r + pd e , i.e., that the 
principal constituents are matter, radiation, and an unknown dark energy, and then infers 
w from the equations of state assigned to each of these constituents. In R h = ct, it is the 
aforementioned symmetries and other constraints from GR that uniquely fix w. 

Both ACDM and R h = ct are Friedmann-Robertson- Walker (FRW) cosmologies, but 
in the latter, Weyl's postulate takes on a more important role than has been considered 
before (Melia & Shevchuk 2012). There is no modification to GR, and the Cosmological 
principle is adopted from the start, just like any other FRW cosmology. However, Weyl's 
postulate adds a very important ingredient. Most workers assume that Weyl's postulate is 
already incorporated into all FRW metrics, but actually it is only partially incorporated. 
Simply stated, Weyl's postulate says that any proper distance R(t) must be the product 
of a universal expansion factor a(t) and an unchanging co-moving radius r, such that 
R(t) = a(t)r. The conventional way of writing an FRW metric adopts this coordinate 
definition, along with the cosmic time t. But what is often overlooked is the fact that the 
gravitational radius, i?h (see Equation 4), which has the same definition as the Schwarzschild 
radius, and actually coincides with the better known Hubble radius, is in fact itself a proper 
distance too (see also Melia & Abdelqader 2009). And when one forces this radius to comply 
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with Weyl's postulate, there is only one possible choice for a(t), i.e., a(t) = (t/t ), where t 
is the current age of the Universe. This also leads to the result that the gravitational radius 
must be receding from us at speed c, which is in fact how the Hubble radius was defined in 
the first place, even before it was recognized as another manifestation of the gravitational 
horizon. 

The principal difference between ACDM and Rh = ct is how they handle p and p. In 
the i?h = ct cosmology, the fact that a(t) oc t requires that the total pressure p be given as 
p = — p/3. The consequence of this is that quantities such as the luminosity distance and 
the redshift dependence of the Hubble constant H, take on very simple, analytical forms 
(as we have already seen in Equation 4). Though we won't need it here, we also mention 
that the evolution of H(z) in the R h = ct Universe goes as 



another very simple and elegant expression that is not available in ACDM. Here, z is the 
redshift, Rh = c/H, and H is the value of the Hubble constant today. These relations are 
clearly very relevant to a proper examination of other cosmological observations, and we are 
in the process of applying them accordingly. For example, we have recently demonstrated 
that the model-independent cosmic chronometer data (see, e.g., Moresco et al. 2012) are 
a better match to R^ = ct (using Eq. 17), than the concordance, best-fit ACDM model. 
Indeed, R^ = ct fits the observed redshift dependence of the Hubble constant, H(z), with a 
reduced Xdof = l-00> compared to the corresponding value xioi = 1-18 for ACDM (see Melia 
2012d). 

In the end, regardless of how ACDM or R h = ct handle p and p, they must both 
account for the same cosmological data. There is growing evidence that, with its empirical 
approach, ACDM can function as a reasonable approximation to R h = ct in some restricted 
redshift ranges, but apparently does poorly in others. For example, in using the ansatz 




(17) 



- 21 - 



P = Pm + Pr + Pde to fit the data, one finds that the ACDM parameters must have quite 
specific values, such as Q m = p m /p c = 0.27 and w^e = — 1, where p c is the critical density 
and Wde is the equation-of-state parameter for dark energy. This is quite telling because 
with these parameters, ACDM then requires Rh(t ) = ct today. That is, the best-fit ACDM 
parameters describe a universal expansion equal to what it would have been with = ct 
all along. Other indicators support the view that using ACDM to fit the data therefore 
produces a cosmology almost (but not entirely) identical to i?h = ct (see Melia 2012c). 

As we shall see below, the results of our analysis of the GRB HD produce very similar 
conclusions to these, i.e., that even though the internal structure of ACDM would appear 
to be quite different from that in = ct (compare Equations 2 and 4), in the end, the best 
fit ACDM model essentially mimics the universal expansion implied by R h = ct. 



In ACDM, the luminosity indicator, and therefore also the distance modulus p, depends 
on the specific choice of parameter values (for Vt m and To directly compare the HD for 
ACDM with that for R h = ct, we will calculate p and <r M using the best-fit model, for which 
Vl m = 0.3 and Q,\ = 0.7. The data and best-fit curve are shown together in the left-hand 
panel of Figure 6. The x 2 f° r this fit is calculated according to 



where p is theoretical value of the distance modulus, and p is measured using the 
Liang-Zhang relation. Also, a- ult is the intrinsic scatter obtained from the joint likelihood 
analysis, and is the error for each realization /i obs of the N data points. Ignoring Hq, 
which does not affect any of these fits, the optimized ACDM model has two remaining 
(principal) parameters, Q m , and so with 28 data points, the reduced x 2 P er degree of 



5. The GRB Hubble Diagram 
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freedom is xL = 46.45/26 = 1.79. 




1 2 3456 1 2 3456 

1+Z 1+Z 



Fig. 6. — Left: Hubble diagram for the GRB sample. The solid curve represents the theory 
(jl in the ACDM model. Right: Same as the left panel, except now the solid curve represents 
the Rh = ct Universe. 

The plot actually gives the impression that the fit is better than this xiof would 
suggest. A closer inspection reveals that 5 data points lie more than ~ la away from the 
best-fit curve. Removal of these data reduces the x 2 considerably, and may be an indication 
that either they are true outliers, or that the errors and intrinsic scatter are greately 
underestimated. 

The Hubble Diagram for the Rh = ct Universe is shown in the right-hand panel of 
Figure 6. Both the data and the best-fit curve were calibrated using the expansion implied 
by this cosmology (see column 5 in Table 2, and Equation 8). A Hubble constant H = 70 
km s _1 Mpc -1 was selected to construct the plot, though it has no bearing on the quality 
of the fit itself. In this case, since we are ignoring Hq in producing the fit, there are 
no remaining free parameters, and the reduced \ 2 P er degree of freedom in R^ = ct is 
x\ $ = 46.79/28 = 1.67. This is a better fit to the data than ACDM, though similar issues 
remain concerning the possible contamination of the GRB sample by outliers and/or the 
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underestimation of errors and intrinsic scatter. To facilitate a direct comparison, these two 
Hubble Diagrams are also shown side by side in Figure 7. 




1 2 3 4 5 6 7 8 



1+Z 

Fig. 7. — A side-by-side comparison of the theoretical curves in ACDM and the Rh = ct 
Universe. 



6. Discussion and Conclusions 

In this paper, we have added some support to the argument that GRBs may eventually 
be used to carry out stringent tests on various cosmological models. Earlier work on this 
proposal had indicated that the spectral and lightcurve features most likely to provide a 
reliable luminosity indicator are the peak energy and a proxy for the jet opening angle, 
which we have taken to be the time at which a break in the light curve is observed. In 
this paper, we have confirmed the notion advanced previously that examining correlations 
among these data can indeed produce a luminosity indicator with sufficient reliability to 
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study the expansion of the Universe. 

A notable result of our work, based on the most up-to-date GRB data, is that a careful 
statistical analysis of these correlations and their optimization points to best-fit parameter 
values in ACDM remarkably close to those associated with the concordance model. We 
have found that the ACDM model most consistent with the GRB Hubble Diagram has 
fl m PS 0.30 and VL\ ~ 0.70. In the concordance model, these values are, respectively, ^ 0.27 
and ps 0.73. 

However, for ACDM the reduced Xdof i s a ^ best approximately 1.79. A close inspection 
of the GRB HD for this model reveals that about 20% of the data points lie at least 2a 
away from the best-fit curve. This may be an indication that some contamination of the 
GRB sample is unavoidable, and that pure luminosity indicators may never be found for 
these sources. Of course, it could also means that we simply have not yet found the ideal 
correlation function, and/or have not yet identified the correct spectral and lightcurve 
features to use for this purpose. On the other hand, it could also mean that we are 
understimating the errors and intrinsic scatter associated with the data. Additional work is 
required in order to better identify the likely resolution to this problem. 

A second principal result of our analysis is that the = ct Universe apparently 
fits the GRB HD better than the optimized ACDM model does. One of our goals with 
this work was to demonstrate the dependence of the data acquisition on the pre-assumed 
cosmological model. This appears to be an unavoidable problem with all cosmological data, 
except perhaps for the cosmic chronometer measurements which are obtained independently 
of any integrated quantity (such as the luminosity distance) that requires pre-knowledge 
of the Universe's expansion history. By calibrating the GRB data separately for ACDM 
and i?h = ct, we have produced a meaningful side-by-side comparison between these two 
cosmologies, showing that the latter fits the GRB HD with a reduced xiof ~ 1-66, noticeably 
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better than the current standard model. Nonetheless, this high value also shows that the 
use of GRBs for cosmological purposes is not yet mature enough to carry out precision 
tests. One may draw the same conclusions here as we discussed earlier for ACDM. 

Another significant result of our study is the remarkable overlap of the two best-fit 
curves in Figure 7. This feature is reminiscent of a similar result from our earlier study 
of Type la SNe, particularly Figure 4 in Melia (2012a). We believe that this is not a 
coincidence because several studies have now shown that ACDM is apparently mimicking 
the expansion history implied by Rh = ct. The most detailed discussion on this issue has 
appeared in Melia (2012c; 2012d). In these papers, we presented several arguments for 
why the optimization of the free parameters in ACDM always seems to indicate an overall 
expansion of the Universe equal to what it would have been in R h = ct. 

Our final comment concerns the implications of this work on the use of Type la SNe to 
study the cosmological expansion at z < 2. There is no question now that any comparative 
analysis between competing cosmologies must be carried with the re-calibration of the data 
for each assumed expansion scenario, particularly when using standard candles that rely on 
integrated quantities, such as the luminosity distance. The Type la supernova luminosity 
cannot be determined independently of the assumed cosmology — it must be evaluated by 
optimizing 4 parameters simultaneously with those in the adopted model. This renders the 
data compliant to the underlying theory. 

Given how much better = ct accounts for the cosmological data, such as the angular 
correlation of the cosmic microwave background (Melia 2012e) and the redshift evolution 
of H(z) (Melia 2012d), not to mention the GRB HD we have studied in this paper, we 
believe it is necessary to produce a Type la supernova Hubble Diagram properly calibrated 
for the Rh = ct cosmology. Only then will it be possible to properly compare the best-fit 
ACDM model directly with R h = ct at z < 2. The payoff from this effort should not be 
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underestimated. We would know for certain whether the Universe is truly now accelerating, 
or whether it continues expanding at a constant rate, as it apparently has been doing from 
the beginning. 
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